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Abstract 

■ 

^•f-} We study the Stephani quantum cosmological model in the presence of a cosmological constant in radiation dominated Universe. In 
the present work the Schutz's variational formalism which recovers the notion of time is applied. This gives rise to Wheeler-DeWitt 

1 'equations which can be cast in the form of Schrodinger equations for the scale factor. We find their eigenvalues and eigenfunctions 
^_by using the Spectral Method. Then we use the eigenfunctions in order to construct wave packets and evaluate the time-dependent 
I expectation value of the scale factor, which is found to oscillate between non-zero finite maximum and minimum values. Since the 
^ expectation value of the scale factor never tends to the singular point, we have an initial indication that this model may not have 

, ^^ singularities at the quantum level. 
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1. Introduction 

In recent years observations show that the expansion 
of the Universe is accelerating in the present epoch [1] 
contrary to Friedmann-Robertson- Walker (FRW) cosmo- 
logical models, with non-relativistic matter and radiation. 
Some different physical scenarios using exotic form of mat- 
ter have been suggested to resolve this problem [2,3,4,5,6,7] . 
In fact the presence of exotic matter is not necessary to 
drive an accelerated expansion. Instead we can relax the as- 
sumption of the homogeneity of space, leaving the isotropy 
with respect to one point. The most general class of non- 
static, perfect fluid solutions of Einsteins equations that 
are conformally flat is known as the "Stephani Universe" 
[8,9] . This model can be embedded in a five-dimensional flat 
pseudo-Euclidean space, which is not expansion-free and 
has non-vanishing density [10,8,11]. In general, it has no 
symmetry at all, although its three dimensional spatial sec- 
tions are homogeneous and isotropic [12]. The spherically 
symmetric Stephani Universes and some of their subcases 
have been examined in numerous papers [9] . So it may be 
important to study the quantum behavior of this model. 

The notion of time can be recovered in some cases of 
quantum cosmology, for example when gravity is coupled 
to a perfect fluid [13,14,15]. This kind of systems are of- 
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ten studied as follows [7,14,15]. First one uses the Schutz's 
formalism for the description of the perfect fluid [16,17], 
second one selects the dynamical variable of perfect fluid 
as the reference time. Finally, one uses canonical quantiza- 
tion to obtain the Wheeler-Dewitt (WD) equation in min- 
isuperspace, which is a Schrodinger-like equation [13]. Af- 
ter solving the equation, one can construct wave packets 
from the resulting modes. The wave packets can be used to 
compute the time-dependent behavior of the scale factor. If 
the selected time variable results in a close correspondence 
between the expectation value of the scale factor and the 
classical prediction (prediction of General Relativity) for 
long enough time, the selected time variable can be con- 
sidered as acceptable. This approach has been extensively 
employed in the literature, indicating in general the sup- 
pression of the initial singularity [7,13,18,19,14,20,15]. 

In the present paper, we use the formalism of quantum 
cosmology in order to quantize the Stephani cosmological 
model in Schutz's formalism [16,17] and find WD equation 
in minisuperspace. In the Schutz's variational formalism 
the wave function depends on the scale factor a, and on the 
canonical variable associated to the fluid, which plays the 
role of time T . Here, we describe matter as a perfect fluid 
matter. Although, this is essentially semiclassical from the 
start, it has the advantage of defining a variable, connected 
with the matter degrees of freedom, which can naturally 
be identified with time and leads to a well-defined Hilbert 
space structure. Moreover, after the universe reaches the 



dust dominated matter, the evolution towards an expo- 
nentially expanding epoch involves a quantum mechanical 
transition associated with some vestige component of the 
original wave function of the Universe. In fact, the classical 
Universe on large scales is based on a quantum mechanical 
background. Particularly, a rapidly oscillating state with 
small amplitude, which is the cosmological influence of a 
wave function vestige, would emerge in the dust dominated 
epoch [26]. 

2. The Model 

The action for gravity plus perfect fluid in Schutz's for- 
malism is written as 

S=\j d^xj—g (R - 2A) + 2 J SxJhh ah K ah 

M dM 

+ J d A x^—g V , (1) 

M 

where K ab is the extrinsic curvature, A is the cosmologi- 
cal constant, and h a b is the induced metric over the three- 
dimensional spatial hypersurface, which is the boundary 
dM of the four dimensional manifold M in units where 
8ttG = 1 [27]. The last term of (1) represents the matter 
contribution to the total action, and p is the pressure. In 
Schutz's formalism [16,17] the fluid's four-velocity is ex- 
pressed in terms of five potentials e, (, £, 6 and S: 



u v = -(e ; „ + C£,„ + 0S iV ), 
I 1 

where p is the specific enthalpy, the variable S is the specific 
entropy, while the potentials Q and £ are connected with 
rotation and are absent in models of the FRW type. The 
variables e and 9 have no clear physical meaning. The four- 
velocity is subject to the normalization condition 



u u. 



= -1. 



(3) 



The metric in spherically symmetric Stcphani Universe 
[10,28,8,12,9,29] has the following form, 



ds 2 



F(t) 



a(t) d (V(r,t) 



V(r,t) dt V a(r,t) 



dt 2 



+ 



a 2 (t) 



(dr 2 +r 2 dn 2 ) 



V 2 (r,t) 

where the functions V(r, t) and F(t) are defined as 
1 



V(r,t) 
F(t) 



1 



-k{t)i 



a(t) 



y/C*(t)a?(t) - k(t) 



(4) 

(5) 
(6) 



Using the line element (4) and the Einstein's equation, one 
can easily show the functions C(t), k(t) and a(t) are not all 
independent, but are related to each other with the follow- 
ing expressions 



Pit) 
p(r,t) 



3C 2 (i) + A 



1 



8ttG ' 
\2C(t)C(t) 



8ttG l 
-3C 2 (t)-A], 



V(r,t)/a(t) 
(V(r,t)/a(t)) 



(7) 



(8) 



where an overdot denotes a derivative with respect to t. 
Note that in the spherically symmetric Stephani models 
and the given coordinate system, the energy density p(t) 
is uniform, while the pressure p(r, t) is not and depends 
on the distance from the symmetry center placed at r = 
0. This is the reason why in such models the barotropic 
equation of state (i.e. of the form p = p{p)) does not exist. 
If, however, we assume some relations between p(t) and 
p(r, t) , this could allow us to eliminate one of the unknown 
functions, e.g. C(t). Hence we are left with two unknown 
functions k(t) and a(t). The first one k(t) plays the role of 
a spatial curvature index, while the second one a(t) is the 
Stephani version of the FRW scale factor. 

Now we consider an observer placed at the symmetry 
center of the spherically symmetric Stcphani Universe. All 
of our physical assumptions will concern the neighborhood 
r w 0. First of all we assume that locally, matter filling 
up the Universe fulfils a barotropic equation of state of the 
standard form 



p(rwO,i) = ap(i). 



(9) 



By substituting the Stephani metric (4) in the action (1) 
and choosing a curvature function k(t) in the form [30] 



(2) k(t)=0a*(t), 



(10) 



and after some thermodynamical considerations [13], the 
final reduced effective action near r w 0, takes the form 



dt 



3^ _ AiV a 3 + 3 pNa 1+ ^ 



,(H) 



+ AT-V« a 3 ^——(e + eS)^ a+1 expf - -) 

(a + l)V«+i v ) \ a) 

where N — F(t)a(t)^ (^(tj)- The reduced action may be 
further simplified by canonical methods [13] to the super- 
Hamiltonian 



12a 



Aa 3 - 3/3a 1+7 + 



p« +1 e s 



(12) 



where p a = —6da/N andp £ = —p u°Na 3 , p being the rest 
mass density of the fluid. Using the canonical transforma- 
tions 



e = e - (a + 1) — , 

Pe 



PT =p a + 1 e s 1 

Pe =Pe, 



(13) 



which are the generalization of the ones used in Ref. [13], 
the super-Hamiltonian takes the form 



3 a 



(14) 



where the momentum px is the only remaining canonical 
variable associated with matter and appears linearly in the 
super-Hamiltonian. 

The classical dynamics is governed by the Hamilton equa- 
tions, derived from Eq. (14) and Poisson brackets, namely 



a = {a,NH} = 
Pa = {p a ,NH} 
-WAa 2 H 



N Pa 
6a 

N 
3Na 



pl + 3N(l + j)f3a-> 



a l+3a 



PT ■ 



(15) 



-3a 



f= {T,NH} = Na 



PT = {p T ,NH} = 0. 



We also have the constraint equation H = 0. Choosing the 
gauge N — a 3a , we have the following equations for the 
system 



T = t, 

3 = (3a-i)^ -^(l+^a 6 ^- 1 



3d 



= -^ T -3/3^+ 1 +Aa 3 + ^. 
n 6a— l a 



(16) 

(17) 
(18) 



Note that the classical equations for the case 7 = +1, in 
Ref. [31], correspond with choosing A = and N = 1. In 
this case (A = 0, N = 1) the constraint equation H = 
reduces to 



iaa z - 3(3a z + a" d > T = 0, 



or 



PT 



(19) 



(20) 



Imposing the standard quantization conditions on the 
canonical variables (p a — > —iJ^,pT — > — *^j0 an< ^ de- 
manding that the super-Hamiltonian operator (14) anni- 
hilate the wave function, we are led to the following WD 
equation in minisuperspace (fi, = 1) 

^-4 + (12Aa 4 - 36/3a 2+7 )* - iUa 1 - 3 "^— = 0. (21) 
oa z at 

According to the equation (16) T = t, can be associated 
with the time coordinate [32,33]. Equation (21) takes the 
form of a Schrodinger equation id^l jdt = . As discussed 
in [19,33], in order for the Hamiltonian operator H to be 



self-adjoint the inner product of any two wave functions $ 
and * must take the form 



(*,*) 



- 3a <5>*^da. 



(22) 



Moreover, the wave functions should satisfy the restrictive 
boundary conditions which the simplest ones are [33,34] 



#(0,i) = or 



dV(a,t) 



da 



= 0. 



(23) 



The WD equation (21) can be solved by separation of vari- 
ables as follows, 



*(<*,*) 



AEt 



V>(a), 



(24) 



where the scale factor dependent part of the wave function 
ip(a) satisfies 

- V" + (36/3a 2+7 - 12Aa 4 )?A - 12Ea 1 - 3a - 



(25) 



and the prime denotes derivative with respect to a. The 
interesting feature of the Stephani model is that the spatial 
curvature is time-dependent. The recent observational data 
shows that our Universe is spatially flat. Moreover, negative 
powers in equation (10) lead to the spatially flat Universe 
in the present epoch. 

We construct a general solution to the WD equation (21) 
by taking linear combinations of the tl> n (a)'s, 



#(a,i) = ^C„V«(a)e 



iE n t 



(26) 



where the coefficients C n (E n ) will be fixed later by choosing 
appropriate boundary conditions, we can compute the ex- 
pected value for the scale factor a for arbitrary wave pack- 
ets, following the many worlds interpretation of quantum 
mechanics [35,36] . We may write the expected value for the 
scale factor a, with regards to the inner product (22) as 



S~a*-s°\*{a,t)?da 
f oo a 1 - 3 «|^'(a,^)| 2 c^a■ 



(27) 



3. The Spectral Method 



To solve the resulting WD equation (25) we use Spectral 
Method (SM) [37,38,39,40] for finding the bound states of 
(25). SM is simple, fast, accurate, robust and stable. This 
method consists of first, choosing a finite domain for the 
approximate solution denoted by 2L, and second taking the 
solution as a finite superposition of the Fourier basis func- 
tions in this domain which satisfy the periodic boundary 
condition. By substituting the expansion into the differen- 
tial equation a matrix equation is obtained. By minimizing 
the first eigenvalue of the resulting matrix with respect to 
the value of spatial domain 2L, the optimized basis func- 
tions can be found. The accurate energy eigenvalues corre- 
spond to the eigenvalues of the resulting matrix with opti- 
mized basis functions. We only examine the bound states 



of this problem, i.e. the states which are the square inte- 
grable. The general equation that we want to solve can be 
written in the form 



E D "i,i U i (~J~) = E A ™,ih{x)Ui (~J~~) 



(38) 



d 2 il){x) -.dip(x) 

^ x )—[^ + M x )— — + M X M X ) 



dx 



ei/>(x), (28) 



Any complete orthonormal set can be used for the SM. We 
use the Fourier series basis. That is, since we need to choose 
a finite subspace of a countably infinite basis, we restrict 
ourselves to the finite region —L<x<L. The value of 
L is determined by requiring that the sought-after eigen- 
functions have compact support in this domain subject to 
the aforementioned optimization. This means that we can 
expand the solution as 



By multiplying both sides of the above equations by 
uv ( m £ x ) and integrating over the x-space and using the 
orthonormality condition of the basis functions, one finds 

L 

•sr-^ /"mir\ 2 f /"mirx\ ... / m'lr. 

Bm,i = ~ 2_, A ™',V [ — ) JUt {—j^- ) h{ x ) u i' [~~[~ 



tp{x) = 
where 

Ul ( 

u 2 ( 



/mirx\ 



2 oo 
i—l m—0 



Or, 



(29) 



L, 



dx 



m TTX 



dx, 



(40) 



. L J 
"miTX\ 



V LR n 
1 



L ) ^/LR 

n 



L 

/mirx\ 



Rr, 



2, m=0, 

(30) 

1, otherwise. D n 



^ y . A m / ,i' 



That is we assume periodic boundary condition. We can 
also make the following expansions 



L 

Ui 

■/ A 



frairx\ , , . / m ux \ 



hi?) 



h(x) 



d 2 xp(x) 
dx 2 

dip(x) 
dx 



^mnx^ 



E c - 



fmTTX\ 



Ui 



\ L 



f 3 (x)ip(x) = ^D m ,i Ui 



(31) 



(32) 



(33) 



= E d ™, 

Therefore we can rewrite Eq. (35) as 



(41) 



E 



+ d n 



eA mji .(42) 



By selecting a finite subset of the basis functions, e.g. choos- 
ing the first 2N which could be accomplished by letting the 
index m run from 1 to TV in the summations, equation (42) 
can be written as 



where B 



m,2 

Cm.i and D r . 



t i are coefficients that can be de- 
termined once fi(x), faix) and fz{x) are specified. By sub- 
stituting and using the differential equation of the Fourier 
basis we obtain 



D A = sA, 



(43) 



n,i 



Bn 



D r . 



>ni —r~ ) ^2^ l A m . l u i\—j-^ 



*V L 



Because of the linear independence of u^ 2 ^), every term 
in the summation must satisfy 



£ An 



(35) 



It only remains to determine the matrices B, C and D. 
Using Eq. (29) and Eqs. (31,32,33) we have 

E /TO7TX\ X /m7T\ 2 (miTX\ , . 

B m ,tu t [—j^- ) = ~2_^ Am >i \~ ) Jn x ) u \—j-^- ) A 36 ) 



mirx\ \ ^ . . d (m-KX 

—) =E^»,*/»(*)^«i(— 



), (37) 



where D is a square matrix with (2N) x (2N) elements. Its 
elements can be obtained from Eq. (42). The eigenvalues 
and cigenfunctions of the Schrodingcr equation are approx- 
imately equal to the corresponding quantities of the matrix 
D. That is the solution to this matrix equation simultane- 
ously yields 27V sought after eigenstates and eigenvalues. 

4. Results 

Since the Hamiltonian commutes with Parity operator, 
eigenstates divide into odd and even categories. We choose 
the odd solutions which satisfy the first boundary condi- 
tion (23). We are free to choose the parameters /3, 7, and 
A. Unlike FRW models in which the bound states merely 
correspond to negative cosmological constants [34] , in this 
model we can find the bound states with negative or posi- 
tive cosmological constant by choosing suitable value of 7. 
In radiation regime (a = 1/3), comparing equations (25) 
and (28) we have 



h(x) = -i 



h{x) = o, 



(44) 
(45) 



/ 3 (x) = 36/?a 7+2 - 12Aa 4 , e = 12E. 

Here, we restrict ourselves to two cases: (/? = 1, 7 = 4, A = 
1) and ((3 — 1,7 = —2, A = —1). In radiation dominated 
regime the expectation value of the scale factor can be writ- 
ten as (27) 



(a)t = 



J °° a\y{a,t)\ 2 da 
f~ma,t)\*da- 



(46) 



Table 1 shows the first 20 odd energy eigenvalues of (25) for 
the two mentioned categories with 10 significant digits. We 
can now construct the wave packets (26) by superimposing 
the resulting eigenfunctions (26). Here, we choose first 20 
eigenfunctions (m — 20) and, for simplicity according to 
Ref. [34], select the coefficients equal to one (C„ = 1) to 
incorporate equally all energy levels. As can be seen from 
the classical equations of motion (17,18), the model has 
singularities at the classical level. At the quantum level the 
probability density of finding the scale factor is (27) 



P(a,t) 



l-3ct 



(47) 



Since the probability density of finding the scale factor at 
a = in radiation regime is zero for odd solutions at all 
times (lim a ^ \^( a ,t)\ 2 = 0), we have a initial indication 
that these models may not have singularities at the quan- 
tum level. Figures 1 and 2 show the behavior of expectation 
value of the scale factor for two mentioned cases in com- 
parison with the classical behavior. Although, in classical 
case the scale factors reach the zero axes, the expectation 
value of the scale factors never tend to the singular point. 
This is depicted in Figs. (3,4) which show the long time be- 
havior of the scale factors. This means that the big bang 
and big crunch phenomena are absent at quantum level. 
Similar properties have been discussed in [13,34,7] for FRW 
cosmological models. 
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Fig. 1. Classical behavior of the scale factor (dashed line) and the 
quantum mechanical expectation value of the scalar factor (solid 
line) for f3 = 1, A = 1, and 7 = 4 in radiation regime. 
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Fig. 2. Classical behavior of the scale factor (dashed line) and the 
quantum mechanical expectation value of the scalar factor (solid 
line) for f3 = 1, A = —1, and 7 = —2 in radiation regime. 
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5. Conclusion 
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In this work we have investigated perfect fluid Stephani 
quantum cosmological model in the presence of cosmolog- 
ical constant. The use of Schutz's formalism allows us to 
obtain a Schrodingcr-like WD equation in which the only 
remaining matter degree of freedom plays the role of time. 
We have obtained eigenfunctions and therefore acceptable 
wave packets have been constructed by appropriate linear 
combination of these eigenfunctions. The time evolution 
of the expectation value of the scale factor has been de- 
termined in the spirit of the many-worlds interpretation 
of quantum cosmology. We have shown that, contrary to 
the classical case, the expectation values of the scale factor 
avoid singularity in the quantum case. 



Fig. 3. The expectation value of the scalar factor for f3 = 1, A : 
and 7 = — 2 in radiation regime for a long period of time. 




Fig. 4. The expectation value of the scalar factor for f3 = 1, A = 1, 
and 7 = 4 in radiation regime for a long period of time. 





/? = 1,7 = 4,A = 1 


P = 1,7 = -2,A = -1 


Ei 


0.7404299830 


3.724923306 


E 2 


2.716898545 


5.221651006 


E 3 


5.467397344 


7.051977995 


E 4 


8.816858136 


9.123953661 


E 5 


12.67642129 


11.39097635 


Ee 


16.98807771 


13.82407714 


E 7 


21.71008635 


16.40315278 


E 8 


26.81052795 


19.11327742 


E 9 


32.26395981 


21.94284703 


Bio 


38.04947945 


24.88253248 


En 


44.14951426 


27.92463946 


E12 


50.54902037 


31.06269279 


En 


57.23492813 


34.29115371 


El4 


64.19574440 


37.60522022 


Eis 


71.42125871 


41.00068165 


Eie 


78.90232098 


44.47381007 


En 


86.63066966 


48.02127735 


Eis 


94.59879650 


51.64009065 


E19 


102.7998386 


55.32754141 


E 20 


111.2274909 


59.08116446 



Table 1 

The lowest calculated odd energy levels for two cases in radiation 
dominated Universe. 
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